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Abstract: We perform nonperturbative studies of the dimensionally reduced 5d Yang- 
Mills-Chern-Simons model, in which a four-dimensional fuzzy manifold, "fuzzy S 4 ", is 
known to exist as a classical solution. Although the action is unbounded from below, a 
well-defined vacuum, which stabilizes at large N, exists when the coefficient of the Chern- 
Simons term is sufficiently small. However, this vacuum corresponds to the "Yang-Mills 
phase", in which the system behaves similarly to the pure Yang-Mills model. In Monte 
Carlo simulations we find that the fuzzy S 4 prepared as an initial configuration either 
diverges or falls into the Yang-Mills phase. Thus the model does not have a "fuzzy S 4 
phase" in contrast to our previous results on the dimensionally reduced 3d Yang-Mills- 
Chern-Simons model, in which the "fuzzy S 2 " can be realized dynamically. 
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1. Introduction 

Fuzzy spheres [1], which are simple compact noncommutative manifolds, have been dis- 
cussed extensively in the literature. One of the motivations comes from the general expec- 
tation that noncommutative geometry provides a crucial link to string theory and quantum 
gravity Indeed Yang-Mills theories on noncommutative geometry appear in a certain low- 
energy limit of string theory [2]. There is also an independent observation that the space- 
time uncertainty relation, which is naturally realized by noncommutative geometry, can be 
derived from some general assumptions on the underlying theory of quantum gravity [3]. 
Another motivation is to use fuzzy spheres as a regularization scheme alternative to the 
lattice regularization [4]. Unlike the lattice, fuzzy spheres preserve the continuous symme- 
tries of the space-time considered, and hence it is expected that the situation concerning 
chiral symmetry [5-19] and supersymmetry might be ameliorated. 

Fuzzy spheres appear as classical solutions in matrix models with a Chern-Simons-like 
term [20-23] and their dynamical properties have been studied in refs. [24-30]. These 
models belong to the class of the so-called dimensionally reduced models (or the large- N 
reduced models), which is widely believed to provide a constructive definition of superstring 
and M theories [31-33]. The space-time is represented by the eigenvalues of the bosonic 
matrices, and in the IIB matrix model [32], in particular, the dynamical generation of four- 
dimensional space-time (in fen-dimensional type IIB superstring theory) has been discussed 
by many authors [34-46]. 

In ref. [28] we have studied the dimensionally reduced 3d Yang-Mills models with the 
cubic Chern-Simons term, which has the fuzzy 2-sphere (S 2 ) as a classical solution. Unlike 
previous works we have performed nonperturbative first-principle studies by Monte Carlo 
simulations. We observed a first order phase transition as we vary the coefficient of the 
Chern-Simons term. For a small coefficient the large N behavior of the model is the same 
as in the pure Yang-Mills model, whereas for a large coefficient a single fuzzy S 2 appears 
dynamically. 
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For obvious reasons it is interesting to extend this work to a matrix model which 
accommodates a /owr-dimensional fuzzy manifold. Among various possibilities, here we 
study the dimensionally reduced 5d Yang-Mills model with the quintic Chern-Simons term 
[21,23], which is known to have the fuzzy 4-sphere (S 4 ) as a classical solution [47]. Various 
aspects of the fuzzy S 4 are discussed in refs. [48-54]. 

The rest of this paper is organized as follows. In section 2 we define the model and 
discuss some properties of the fuzzy S 4 solution. In section 3 we present our results of the 
numerical simulation. Section 4 is devoted to a summary and discussions. 

2. The model and the fuzzy S 4 

The model we study is defined by the action 

S = N (-^tr[A^A u \ 2 -^\e^ paT tr{A p A v A p A a A T )^ , (2.1) 

where A^ (/x = 1, ■ ■ ■ , 5) are N x N traceless hermitian matrices, and e pupaT is a rank-5 
totally-antisymmetric tensor with £12345 = 1. 
The classical equation of motion is given by 

[Ay, [A p} A u ]] + Xe pupiJT A u A p A a A T = , (2-2) 



which is known to have a solution of the form 

1 

2 



A p = \aG^ (2.3) 



where a is given by 



a = 7 rr (2-4) 

(n + 2) A y 1 



with some positive integer n, and G p is defined as 



G 



n 

p = V ( 14 ® • • • ® 14 ®r M <g> i 4 ® • • • ® 14) • (2.5) 

L — ' V s v / v v sym 



sym 

3 j-1 n-j 



Here the tensor product <S> with the suffix "sym" represents the so-called symmetric tensor 
product (See e.g. ref. [55] for details), and T p are 4x4 hermitian matrices which obey the 
5-dimensional (Euclidean) Clifford algebra 

{T p ,T l/ } = 2 5 pu . (2.6) 

An explicit representation of the 5d gamma matrices is 

T a = a a ®a 1 (a = 1,2, 3), (2.7) 

r 4 = i 2 ® o 2 , (2.8) 
r 5 = i 2 ® a 3 , (2.9) 
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where a a (a = 1, 2, 3) are the Pauli matrices . According to the definition of the symmetric 
tensor product, the size N of the matrix Gp in (2.5) is given by 

N = l( n +l)(n + 2)(n + 3) . (2.10) 

The solution (2.3) can be regarded as the fuzzy S 4 [21,23]. 1 Unlike the fuzzy S 2 case [28], 
the fuzzy S 4 solution exists only for particular values of N given by (2.10). 

From the algebraic point of view, the fuzzy S solution (2.3) differs from the fuzzy 
S 2 in that Gp is not closed with respect to the commutator. Let us therefore define the 
commutator as Gp V = [Gp, G v \. Then Gp and Gp U satisfy the following relations [47,48,55] 

GpGp = n(n + 4)l N , (2.11) 

G^Guv = -16 n(n + 4) 1 N , (2.12) 

[GV, G p ] = {-SupGv + S up G^) , (2.13) 

[Gp,u, G px ] = (S I/p Gp X + bp X G vp — bp P G vx — b vx Gp P ) . (2-14) 

From eq. (2.11) we obtain 

^(^) 2 = i? 2 l^, (2.15) 

R 2 = i a 2 n(n + 4) , (2.16) 

which means that R represents the "radius" of the fuzzy S 4 . Note that the radius R is 
inversely proportional to A, the coefficient of the Chern-Simons term, unlike the fuzzy S 2 
case [28]. The matrix Gp also satisfies the self-duality condition 

e^ P arG u G p G a G T = 8(n + 2)Gp. (2.17) 

Using these formulae, it is easy to check that (2.3) is a solution to the classical equation of 
motion (2.2). The value of the action for the classical solution (2.3) is given as 

S d = ^n(n + 4)a 4 N 2 , (2.18) 

which increases with a and hence with the radius R of the fuzzy S 4 . This is again in 
contrast to the fuzzy S 2 case [28], where the corresponding action has an overall minus 
sign, and therefore it decreases with the radius of the fuzzy S 2 . 

In order to reveal some crucial properties of the present model, let us consider a one- 
parameter family of configurations 

Ap = ^Gp, (2.19) 



1 More precisely, the solution (2.3) corresponds to the homogeneous space SO(5)/U(2), which is locally 
S 4 xS 2 . Here the fuzzy S 2 is naturally considered as an internal space attached to each point on the fuzzy 
S 4 [21,23]. 
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where £ = a corresponds to the classical solution (2.3). The corresponding value of the 
action (2.1) for each £ is given by 

S(0 = N 2 n(n + A) - J-£ 5 ) . (2.20) 

We find that the action can take arbitrarily negative value by increasing £. Thus the 
classical action is not bounded from below and the path integral over obviously diverges. 
Moreover, the point £ = a actually gives a local maximum of S(£). Thus the solution (2.3) 
is merely a saddle point of the classical action. These properties are in striking contrast 
to the situation in the fuzzy S 2 case [28], where the single fuzzy S 2 is a minimum of the 
action (at least locally), and the path integral converges for TV > 4 [56]. These differences 
come from the fact that the Chern-Simons term in the present model is quintic and it 
has a higher power in than the Yang-Mills term, whereas in the previous model the 
Chern-Simons term is cubic and it has a lower power in A^ than the Yang-Mills term. 

Even if the path integral diverges for A ^ at finite N, there is a possibility that the 
model possesses a well-defined vacuum, which stabilizes at large N. A well-known example 
is the one-matrix model 

5^ = ivQtr0 2 -i 5 tr^ , (2.21) 

where eft is an N x N hermitian matrix. Although the action is unbounded from below, 
this model has a sensible large- N solution [57] for g 2 < g 2 = j^g - Similarly, our model 
(2.1) does possess a well-defined vacuum for sufficiently small A, as we will see in the next 
section. We will also examine the possibility that the fuzzy S 4 , which is unstable classically, 
may get stabilized due to some quantum effects. 

3. Monte Carlo simulations 

We address the issues raised in the previous section by Monte Carlo simulation. For an 
updating procedure, we apply the heat bath algorithm developed in ref. [58] (See the 
appendix A of ref. [28] for more details). We take the fuzzy S 4 solution (2.3) as the 
initial configuration, and carry out simulations for various A at N = 4, 10, 20, 35, 56, which 
corresponds to n = 1,2,3,4,5, respectively, due to (2.10). Since the qualitative behavior 
turns out to be the same for any N, we will show explicit results only for N = 56. In 
order to probe the geometrical structure of the configurations, we consider the eigenvalues 
Xj (j = 1, • • • ,N) of the quadratic Casimir operator [28] 

Q = £(V • (3-1) 

n=i 

In figure 1 we plot the eigenvalues {Xj} and the normalized action ^ for a = 0.4 
against the Monte Carlo time r (or more precisely, the number of "sweeps" in the heat 
bath algorithm [28]). For the initial fuzzy S 4 configuration (r = 0), all the eigenvalues 
coincide with (2.16) and the action is given by (2.18). After the first sweep, the eigenvalues 
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spread out, obviously for entropical reasons, and correspondingly the action increases. The 
eigenvalues and the action are stable for several sweeps, but change drastically at r = 6. 
At r = 7, which is not included in the figure, the maximum eigenvalue reaches 2.0 x 10 4 
and the action becomes = —1.7 x 10 6 . Subsequently all the eigenvalues blow up and the 
action takes further negative values until the run terminates due to overflow. We observed 
similar behaviors for a < 0.5 (i.e., large A). 




Figure 1: The eigenvalues of the quadratic Casimir operator Q (left) and the normalized 

action (right) are plotted against the Monte Carlo time t for a — 0.4 at N — 56. 

For a > 0.5 (i.e., small A), on the other hand, the system turns out to be well-behaved. 
In figure 2 we plot the eigenvalues {Xj} and the normalized action against the Monte 
Carlo time r. The eigenvalues spread out after the first sweep, and then they all decrease 
keeping their extent almost constant. The action goes up after the first sweep, and then 
goes down as all the eigenvalues decrease. This behavior can be understood qualitatively 
from the action (2.20) for the rescaled configuration (2.19). The action after thermalization 
turns out to be almost independent of a. This is a property of the "Yang-Mills phase" [28], 
in which the effect of the Chern-Simons term is negligible and the model behaves similarly 
to the pure Yang-Mills model (A = 0, i.e., a = oo ) [58]. 

We also calculate the eigenvalue distribution defined by 

1 N 

/(x) = -^-A i )) . (3.2) 

For the fuzzy S 4 solution, it is given within the classical approximation as 

f(x) = 5(x - R 2 ) , (3.3) 

where R 2 defined by (2.16) is proportional to a 2 . Figure 3 shows the results for a = 
0.5, • • • , 1.0. The results are almost independent of a, which continues to be the case for 
even larger a. This is again a property of the "Yang- Mills phase" . Since the empty region 
at x < 1 can be understood from the uncertainty principle, the geometry should rather be 
considered as that of a solid ball [28]. 
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Figure 2: The 56 eigenvalues of the quadratic Casimir operator Q for a = 1.0 (left) and the 
normalized action for a = 0.5, • • • ,1.0 (right) are plotted against the Monte Carlo time r at 
N = 56. 



We have also simulated the model 
using = as the initial configura- 
tion. In this case we are not restricted 
to the values of N given by (2.10), so 
we tried N = 12, 14, 16 as well. For any 
value of N, we find that the system ei- 
ther diverges or falls into the Yang- Mills 
phase. 

4. Summary and discussions 
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Figure 3: The eigenvalue distribution f(x) is plot- 
ted for a = 0.5, • • • , 1.0 at N = 56. 



In this paper we reported on a first at- 
tempt to extend our previous work [28] 
on the fuzzy S 2 to a four-dimensional 

fuzzy manifold. We considered the fuzzy S 4 realized as a classical solution in the dimen- 
sionally reduced 5d Yang-Mills-Chern-Simons model. We find various differences between 
the present model and the previous one in spite of their similarity. In particular we have 
observed that the fuzzy S 4 prepared as the initial configuration either diverges or falls into 
the Yang-Mills phase. Thus we conclude that the present model does not have a phase in 
which the fuzzy S 4 becomes stable. 

This negative result is essentially due to the fact that the quintic Chern-Simons term 
has higher powers in than the Yang-Mills term. 2 In this regard we note that there 
are other models which accommodate four-dimensional fuzzy manifolds such as CP 2 and 
S 2 x S 2 with some kind of a cubic Chern-Simons-like term. We would like to report on 
these cases in the forthcoming publications. 

2 Fuzzy spheres appear as classical solutions also in dimensionally reduced Yang-Mills models with a 
tachyonic mass term [59]. We obtained clear indications from Monte Carlo simulations that such systems 
are ill-defined. 
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